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METHODS AND APPARATUS FOR 
IDENTIFICATION AND IMAGING OF SPECIFIC 
MATERIALS 



BACKGROUND OF THE INVENTION 

10001] This invention relates generally to methods and apparatus 
usef ul for detuning the composition of matenals under study and -P-ady 
to methods and apparatus for analyzing and/or imaging specie matenals m objects 
under study. 

[0002] Because many town CT detection systems do no. provide 
energy teso.ntion, i, is not possMe to provide materia, characterization informatton 

density. As a resnlt, known computed tomograph.c (CT) .mages do no 

latenL ma. have similar density hut different atomic numbera, and .mageamt. -k 

suhstandany uniform even though an objec, under study haa v_ m « 

composition, In addition, beam-hardening artifacts, such as no„-un .formny 

and sneaking can resutt. from the non-linear relationship hetween 

^ path iengms fo, polychromatic x-ray heams in CT imaging system, Als. known 

CT imaging systems do no, provide c,uan,i,a,ive image values. Instead, me same 

material a. different locations can show different CT numbers. 

[0003] At least one known dual energy decomposition algorithm is 
known .ha. represents material-specific characteristics as a two-parameter has.s set 
Is, b y encoding each of rhese parameters individually, two separate unages can he 
form d using a CT system. In a, leas, one known system, a smgle shce unage .s 
Z 1 using a sing,e sfice CT detector system, using two different x-ray beam 

used to obtain scan .he same slice of an object. The two scans are no, perform 
simuhaneously, bu, instead are performed a, shghtty different fitnes, e.g., se q uen„a„y. 
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„ another known sys.cn, energy sensittve scanning is performed by us,„g an ene^y 
live de.ec.or system such as a photon coun.ing de.ec.or. In either case, .he w 
Te gy dependen. da, se.s are used with an appropfiate m,er,a, decompo»..on 
JL I prodr.ee two .nrages, each represent one of .he two basts nra.enais. 

,0004] in wo basis ma.enal decomposition images produced by 
tmaging systems, the imaging vaiue for each pixe, in an 
m a,eriai denstty for me corresponding basis materia.. Any matena, o, 
basis ma.ena,s »H. show up in both images, with the tmage ptxe, bemg 
proportional to dre density of the non-basis ma.enal. 

,0005] Using known two basts material decomposition algorithms, 
my nraterta, Cher than the two basis materials appears in both 
JL an incorrect density. This contamination reduces .he vrstbth.y of h basts 
Ha, in .he images, and also results in denstty errors in q uan.rf,ca.ton 

applications. 

BRIEF DESCRIPTION OF THE INVENTION 

,0006] Some configurations of the present invention therefore 
, ... ■ _ oh ! ect The method includes acquinng 
provide a method for analyzing matenals man object. 

x-ray projection data of .he object a, high energy and a. low energy for a plurahty of 

determine materia! densities a, each pixel for two se.ected bas,s matenals. A 
tuition of an object at each pixel is determined utilizing a dctermmed mappm 
Zteria, density regions for the .wo selected basis m,eHals. An ,mage mdtcattve 
of I composifion of the object is displayed utifizing ,he determined compostbon. 

,0007] Anottter aspect of .he present invention provides a method for 
anting materials in an object .ha, includes acquiring x-ray projection data of the 
S a high energy and a, low energy for a plurafity of views The ^d 
p4ct,on data is ufil.zed in a ma.eria, decomposifion ,o de.erm.ne — ~ 
a, Lh Pixel for two selected basis material, A determtned mapptng of malena, 
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densi.y regions for .he .wo selec.ed basis ma.ena.s is utilize* » fiher pixels of an 
image of .he obiec. con.spond.ng to one o, more compositions of in.eres,, and an 
, mag e ind.cative of .he locations of composition of interes, of .he objee. is d.splayed. 

[0008] Still other aspects of the present invention provide an 
apparatus for analyzing materials in an object The apparatus includes an x-ray source 
aid a de.eo.or configured to acquire projection data a. high and low energ.es for a 
polity of views. The apparams also includes a compu.er, a s.orage dev.ee an a 
display The apparatus is configured .o acquire x-ray projection da.a of the objec. a 
Z -rgy and at low energy for a plurality of views utilizing the x-ray source and 
aTde.ec.or. The ap P ara.us is farmer configured .o process me acqmred x-ray 
projection da.a utilizing the compu.er and .he storage device .o de.«rm,ne ma,ena. 
densities a. each pixel for two se.ected basis materials and to determine a compos. hon 
of an objec. a. each pixe, utilizing .he computer and a defined mappmg of mafenal 
density regions for «he two selected basis materials in .he storage device. The display 
is confignred ,o disp!ay an image ind.ca.ive of .he composition of .he objec. utihzmg 
the determined composition. 

[0009] in yet another aspect of the present invention, there is 
provided an apparatus for analyzing material, in an objec.. The apparatus includes an 
lay source and a de.ec.or configured ,o acquire projection d,a a, h.gh and low 
energies for a plurality of views, a computer, a storage device, and a d.splay. The 
applms is configured <o acquire x-ray projection da. of me objec. a. h.gh en g^ 
Id a, low energy for a plurahty of views ntilizing the x-ray source me sard dc. Cor. 
Xhe apparams is also configured to process me acquired x-ray projection data uhhzmg 
«. compulcr and me storage device to determine materia, densities a. each p.xe, for 
,wo selected basis materials and to ntilize the computer and a determined mappmg of 
ma.eria! density «gions for .he W o selected basis materials in the storage dev.ee to 
fiMer pixels of an image of the objec, copending .o one or more compos.tions o 
ta . Jt. The disp,ay is used to d,sp,ay a, image indicative of tire locations of 
composition of interest of the object. 
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[0010] It will thus be appreciated that configurations of the present 
invention are capable of proving quantitative material information. Some 
configurations of the present invention further prov.de improved target material 
visibility and quantification accuracy and/or provide decomposition into more than 
one material. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] Figure 1 is a pictorial view of a CT imaging system 

embodiment. 

[0012] Figure 2 is a block schematic diagram of the system 
illustrated in Figure 1. 

[0013] Figure 3 is a flow chart of a configuration of the present 
invention used to prepare images from projection measurements at high and low 
energies for a plurality of views of an object. 

[0014] Figure 4 is a flow chart of a method useful for obtaining a 
calibration matrix used in the configuration of the present invention represented by the 
flow chart of Figure 3. 

[001 5] Figure 5 is a flow chart representative of configurations of the 
present invention that produce an image indicative of the composition of an object or 
an image showing locations of one or more composition of interest within the object. 

[0016] Figure 6 is a graph illustrative of a first type of rectangular 
segmentation of material density regions useful in configurations of the present 
invention. 

[0017] Figure 7 is a graph illustrative of a radial segmentation of 
material density regions useful in configurations of the present invention. 
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[0018] Figure 8 is a graph illustrative of a second type of rectangular 
segmentation of material density re gl ons useful in configurations of the present 
invention. 

[0019] Figure 9 is a multiple material decomposition iodine image of 
blocked lumens in a phantom produced by a configuration of the present invention. 

[0020] Figure 10 is a multiple ma.erial decomposition calcium image 
of Hocked lumens in the phantom of Figure 9, also produced by a configuration of the 
present invention. 

[002 1 ] Figure 11 is a conventional, prior art CT image of the blocked 
hm e„s in the phanlom of Figure 9, presented for comparison with Figures 9 and .0. 

DETAILED DESCRIPTION OF THE INVENTION 

[0022] Example embodiments of methods and apparatus for multiple 
materia, decomposttion for energy discnmination are described below in more detatl. 
A technical effect of the methods and apparatus described herein include a. leas, one 
of facilitating the extraction of information concemtng a object utilizing a CT nuagmg 
syst em, generation of images representative of a pmrality of different matenals 
contained within such an object, and automated identification of matenals contamed 
within an object. 

[0023] In some known CT imaging system configurations, an x-ray 
SO urce projects a fan-shaped beam which is coUimated to He within an X-Y plane ,« * 
Cartesian coordinate system and general.y referred to as an '•imaging plane . The x- 
ray beam passes throng an object being imaged, such as a patient. The beam, after 
being attenuated by the object, impinges upon an array of radiation detectors. The 
intensity of the attenuated radiation beam received a. the detector array is dependent 
upon the attenuation of an x-ray beam by the object. Each detector element of the 
array produces a separate electrical si^a, that is a measurement of the beam tntensfty 
at Ore detector location. The intensity measurements ftom all the detectors are 
acquired separately to produce a transmission profile. 
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m cvcr^Ttr? the x-ray source and the 
[0024] In third generation CT systems, me x y 

detector array are rotated with a gantry within the imagmg plane ana around the oh,ec. 
t0 b e imagefi snch tha, the angle a. which ,he x-ray beam mtersects the ob^t 
constantly change, A group or x-ray — n measurements , ie„ P.— 
from the detector array a. one gantry ang,e is referred to as a v,ew A scan of 
ohjec, conrprtses a se, of views made a, different gantry angies, or vtew a„g,es, dunng 
one revolution of the x-ray source and detector. 

,0025] In an axial scan, the projection data is processed to construct 
« image that corresponds to a two-dimenstonal shoe taken through the object One 
leZfor reconstructing an image from a se. of projection data is referred to t„ he 
metnoa tor re nr o C ess converts the attenuation 

art as the filtered backprojecoon techntquc. Thts process con 
m e_,s fiom a scan rnto integers cafied "CT numbers- or "Hounsfi e,d urn. 
(HU) , which are used to contro. the brightness of a correspond ptxel on a cathode 
ray tube display. 

,0026] To reduce fire total scan time, a "helical" scan may be 
performed. To perform a "be.ica," scan, fire pafien, is moved while fine dabs for the 
prescribed number of slices is acquired. Such a system generates a stngle hehx from a 
ZZ helica. scan. The hehx mapped on, by the fan beam yields proton data 
from which images in each prescribed slice may be reconstructed. 

[0027] Reconstruction algorithms for helical scanning typically use 
helical weighing algorithms that weigh, fire collec.ed data as a Junction of view ^ 
and detector channel index. Specifically, prior to a filtered backprojectton process 

the gantry angle and detector angle. The weighted data is then processed to ge era. 

CT numbers and ,o construe, an image mat corresponds to a two-d.mens,onal shoe 

taken through fire object. 

[0028] To further reduce the total acquisition time, multi-slice CT 
„as been introduced, hr mn.ti-s.ice CT, mu.tip.e rows of projection data are ^ 
simult aneous.y a, any time instant. When combined with hehca. scan mode, the 
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system genera.es a single helix of eone beam projection data. Similar to the single 

slice helieal, weighting scheme, a method can be derived ,o multiply the weight wtth 

the projection data prior to the filtered backprojeetion algorithm. 

[0029] As used herein, an element or step recited in the singular and 

proceeded with the word "a" or "an" should be understood as not excluding plural said 
elements or steps, unless such exclusion is explicitly recited. Furthermore, references 
,„ -one embodiment" of the present invention are no, intended to be interpreted as 
excluding the existence of additional embodiments that also incorporate the rectted 
features. 

[0030] Also as used herein, the phrase "reconstructing an image" is 
not intended to exc.ude embodiments of the present invention in which data 
representing an image is generated bn. a viewab.e image ,s no,. However, many 
embodiments generate (or are configured ,o generate) a, leas, one viewable .mage. 

[0031] Referring to Figures 1 and 2, a multi-slice scanning imaging 
system, for example, a Computed Tomography (CT) imaging system 10, is shown as 
including a gantry 12 representative of a "third generation" CT tmaging system. 
Gantry ,2 has an x-ray tube 14 (also called x-ray source ,4 herein) ,ha, projects a 
beam of x-rays 16 toward a detector array 18 on the opposite side of gantry 12. 
Detector array 18 is fotmed by a plurality of detector rows (no, shown) inc.ndmg a 
plurality of detector elements 20 which together sense me projected x-rays ,ha, pass 
through an object, such as a medical patient 22 between array 18 and source 14. Each 
detector element 20 produces an electrical signal Una, represents me intenstty of an 
impinging x-ray beam and hence can be used to estimate me attenuation of ,he beam 
as i, passes through objec, or palien, 22. During a scan to acquire x-ray projection 
data gantry 12 and the components mounted therein rotate about a center of rotation 
24 'Figure 2 shows only a single now of detector elements 20 (i.e., a detector row). 
However, multi-slice detector array 18 includes a plurality of parallel detector rows of 
detector elements 20 snch mat projection data corresponding to a plurality of quast- 
parallel or parallel slices can be acquired simultaneously during a scan. 



-7- 



134431 



Express Mail No. EL 977935786 US 



[0032] Rotation of components on gantry 12 and the operation of x- 
ray source 14 are governed by a control mechanism 26 of CT system 10. Control 
mechanism 26 includes an x-ray controller 28 that provides power and timing signals 
to x-ray source 14 and a gantry motor controller 30 that controls the rotational speed 
and position of components on gantry 12. A data acquisition system (DAS) 32 in 
control mechanism 26 samples analog data from detector elements 20 and converts 
the data to digital signals for subsequent processing. An image reconstruct 34 
receives sampled and digitized x-ray data from DAS 32 and performs high-speed 
image reconstruction. The reconstructed image is applied as an input to a computer 
36, which stores the image in a storage device 38. Image reconstruct 34 can be 
specialized hardware or computer programs executing on computer 36. 

[0033] Computer 36 also receives commands and scanning 
parameters from an operator via console 40 that has a keyboard. An associated 
cathode ray tube display 42 allows the operator to observe the reconstructed image 
and other data from computer 36. The operator supplied commands and parameters 
are used by computer 36 to provide control signals and information to DAS 32, x-ray 
controller 28, and gantry motor controller 30. In addition, computer 36 operates a 
table motor controller 44, which controls a motorized table 46 to position patient 22 in 
gantry 12. Particularly, table 46 moves portions of patient 22 through gantry opening 
48. 

[0034] In one embodiment, computer 36 includes a device 50, for 
example a floppy disk drive, CD-ROM drive, DVD drive, magnetic optical disk 
(MOD) device, or any other digital device including a network connecting device such 
as an Ethernet device for reading instructions and/or data from a computer-readable 
medium 52, such as a floppy disk, a CD-ROM, a DVD or another digital source such 
as a network or the Internet, as well as yet to be developed digital means. In another 
embodiment, computer 36 executes instructions stored in firmware (not shown). 
Computer 36 is programmed to perform functions described herein, and as used 
herein, the term computer is not limited to just those integrated circuits referred to m 
the art as computers, but broadly refers to computers, processors, microcontrollers, 
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circuits, and other programmable circuits, and these terms 

herein. AUho Ug h the specific embodimen, mentioned above refers , a timd 
, „„ CT system the methods described herein equally apply to fourth 

cCems Cslnar, detector and x-ray source, Add— ,« , — * 
, hat the benefit of the invention accrue to imaging modahfes other than CT 
tlionaUy although the herein described methods and apparatus are described rn a 
^1 1 X it ifeontemplated that the benefus of the invention accrue to non 

*— - - - *— tt« 

setting or a ttansporiation setting, sueb as, for example, but no. hmt.ed to, a baggag 
scanning system for an airport or other transportation center. 

,00351 To provide material decomposition, some configurations of 
CT detector anay include two differen, x-ray beam filters (no, shown), or x-ray tube 
I temted Irwo differen, x-ray tube vouages <xV,s>. * «- «- « 
vohages are used ,o ob,ain scans of ,he same slice or shccs of ob3ec, 22 The two 
™ no, perforated simuttaneousiy in tirese configurations. Sequenna, scans 

slice or slic es are satisfactory in studies in which littte movemeu, ,s expec.ed 
bleen scans. Some omcr configurations prov.de an energy sensitive 
iZTly .S, such as a pho.on counting detec.or, for energy sensitive scaumng 
Conflations a„o„ the simuhaneous acquisition of scans a, differen, energy 
Tier Z the two energy depcnden, data sc. are used with an appropn* 

two basis materials. 

,0036] In some configurations of the present invention, a CT imaging 
system 10 acqutrcs dua, energy projection measurement da,a aud processes ,he dara ,o 
system , q . Ba ,h n f these images represents a different material 

provide two or more images. Each of fine ,mag P 
signature of object 22 under study, energy drscnmmation CT is used 
Igt system is used ,0 separa.e.y detect tivo regions of pbo.on energy specp 
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r 0 037] More particularly, in various configurations of the present 
Mention, the measured x-ray attenuation for a material at two or more x-ray energ.es 
rimpos^ 

Z 7e) is a measured —on coefflcien, to, ma t er,a, m (an « materia,) a, 
energy E ; 

c is a linear coefficient for a first material A; 

i m is a linear eoeffieient for a second material B; 

„" (£) is an attenuation eoeffieient for material A at energy E ; and 

„',(£) is an attenuation eoeffieient for material B a. energy E . 

,0038] Thus, in some configurations of the present invention and 
referring to flow char, 100 of Figure 3, a technical effect of imaging system 10 ts 
"d hy a user by operating imaging system 1 0 to obtain predion measurements 
TZ l >ow en rgies for muUip.e views a, ,02. Computer 36 and/or srorage 
3 7l used in Le configurations ,o generate materia, decomposition .mages 
72 or ,wo se.ee.ed materiaU using any suitahie method, such as one ustng tab 
W Tnps hr such confiscations, a two-dimensiona, table is stored m memory 
ZZ« eompnrer 36. This two-din— rs referenced usnrg rivo ,^e 
one of which corresponds to pixe, vahres of the ,ow energy rmage and flte other 
Ich co^sponds ,o va,ues of corresponding ,xe,s of the high energy .mage. Th 

„ corresponding to densities of materia.s A and B. The values p., and p., are 
« as densities of A and B, respective,* a, the pixe, .ocation. (The arrangement 
Z ,ah,e in memory is ,eft as a design choice. For examp,e, some configwafl 

p., , rather man a single table in memory tha, provides both entries.) 
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[0039] In some medieal imaging applications, water and iodine can 
usefully be selected as ma.ena, A and meter,, B, respective^, bu, neither the 
, n ven,io„ i.se, f nor its medical imaging applicafions are l.mi.ed to tb.s selecon of 

inferred densities of each malerial are reconstructed in some configurations a. ,04. 

[0040] In some configurations and referring to flowchart 200 of 
Hgure 4, the contents of the table in memory refen-ed to above is predetermined by a 
user, a suppher, or anolher entily by opening an Imaging system 10 or ofcer x my 

eLgy data as a function of thickness of matena, A and of ma,ena, B a. 202. Th 
virion ,s used to generate a calibration matrix a, 204, and inverted a, 206 o 
generate ,he maleria, decomposition .ookup rab,e. Alflrough a plummy o 
Lasurements are reauired a, 202 to obtain me needed data, the cahbrauon procedure 
representor, by flowchart 200 need on,y be performed once for malenals A -M 
J„ be noted ma, eonf.gura.ions of the present invenlion may be provided wt* 
programmed lab.es, or ,ab,es Iha, may be programmed from machme readable 
me ia, such as media 52. to some configurations, the lab.es may be determined by a 
user of imaging system .0 a, a ,ime chosen by the user, upon me usef s recogmfon of 
the need for a particular table in a particular s.tuanon. 

[0041] In various configurations of the present invention and 
referring to flowchart 300 of Figure 5, .he leehnical effect of imaging system 10 ,s 
achieved by a user operaung imaging sys.em .0 to acuuire projection measurements 
or a plurality of views of an object 22 a, high and low energies. At 304, irrespeCve 
of me ac.ua, composition of ohjec. 22, a material decomposition is performed 
m e .able derived for materia* A and B. The two .inear coefficenls derived from flu 
Huear combination are dislinclive for me various maleria.s and will characterize*, 
effective atomic number of a materia,. Using the two linear coefficients *us denv«l 
e„ab,es configurations 10 of the present inven.ion .o reconstruc. a materra. de„,.y for 
anv targe, 22 materia, and ,o encode an image. This reconstruction creates separa 
cross-seciona, or vo.umelric da,a sea for each ,arge, 22 ma,eria,. As a resu,,, 
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visibility and quantification accuracy of composite and spatially vanable objects is 
improved. 

[0042] More particularly, in CT system 10, a CT number in a 
reconstructed image represents a linear attenuation coefficient: 

where p is a materia, density of objeet 22 having nnits g/cm', „. is a mass 
attenuation coefficient having units cm'/g, and a is a linear attenuation coefficen, 
having units em '. Thus, if the materia, of objeet 22 is known, one can detenrune the 
density p. 

[0043] Ignoring k-edge effects, x-ray attenuation for any matenal is 
decomposed into a linear combination of the attenuation for two basis materia* A and 
B These coefficients are different for materials that have different effecttve atomre 
number. Using these different coefficients, conflgurations of the present inventton 
utihzing muttiple material decomposition to reconstruct tbe matenal density .mage for 
any target materials. Thus, 

where p t is a materia, density of basis materia. A having units g/cm', p., ts a 

nraterial density of basis materia. B having units g/cm\ „, is a mass attenuation 

eoeffieien, of material A, and p, is a mass attenuation coefficient of materia, B. 

Referring to Figure 3, PmA and can be represented in a radial coordinate system 

/ 2 , n 2 an( i an angle written 
by a radial distance wntten r m - V Am + P«* 

6 - tan-<p / A.) • As noted above, combinations of PmA and p mB derived from a 
serial decomposition are distinctive for various materials including materials other 
than A and B, as are r„ and 9 . 

[0044] Therefore, using a predetermined mapping of material density 
regions, a composition of an object a, each pixel (or BUer pixels corresponding to one 
or more compositions of interest) is detennmed by image reconstntctor and/or 
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computer 36 a, 306. The determined composition information is then used, m some 
configurations, to prepare an image indicative of the compositton of ,he ob.ec, (or of 
the locations of compositions of interest in the object) at 308. 

[0045] More particularly and referring to Figures 6, 7 and 8, multiple 
materia, decomposition is achieve, by firs, acquinng dua, energy measurement 
projection data. Any known two materia, deeompositton method is then applied ,o the 
dual energy measurement projection data. The parameters p. 4 and p., serve as x 
and y coordinate values in a Cartesian coordinate system, or equivalent!* r. and 0 
serve as coordinate values in a radial coordinate system. In Figures 6, 7, and 8, for 
example, water is basis materia! A and iodine is basis materia! B, so that p. . an x- 
axis coordinate and p, is a y-axis coordinate. In this coordinate system, a 
decomposition of a pure Ca (calcium) object has coordinates x-1.188, y=0.1088. 
Mixtures of ealeium in water fall on a line from the pure calcium coordinates to x-1, 
v=0 depending upon the concentration of Ca. Mtxtures of an and ealeium fall on a 
fine from the pure ealeium coordinates to x=0, y=0, depending upon the relative 
amounts of calcium and air a, a given pixel point. A water-iodine mixture has 
coordinates on a line from x-0, y-4.94 to x=l , y-0, depending upon concentrafon. 

[0046] To take advantage of the loci of various materials in the x-y 
coordinate system, regions of the plane defined by the x-y coordinates are segmented 
in some configurations of the present invention and assigned ,o a material that ,s moa 
Hk* represented by mat segment. For example, in Figure 6, calcium, water, an, and 
iodine are materials of interest in a medical image, and rectangular segments SI S2 
S3 S4 S5 S6 S7.S8, and S9 are identified. Segments SI, S2, and S3 are considered 
"out of range" and are no. associated with materials likely to be in the object being 
imaged Segment S4, is indicative of large iodine concentration, although an .mage 
pixel is unlikely to fall into this segment. Segment S5 is indicative of a calctum-atr 
mixture. Segment S6 is indicative of a high concentrafon of ealeium. Segment 87 
does no, definitively indicate a composition (a pixel falling into mis segment .could, 
for example, be indicative of mostly air with some ealeium, or a relatively htgh 
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concentration of iodine in water). Segment S8 also does not definitively indicate a 
composition. Segment S9, however, is indicative of water, with either a low 
concentration of 10 dine or no more than a moderate concentration of calcium, so a 
pixel falling into Segment S9 is interpreted as being indicative of water. 

[0047] Some configurations of the present invention use radial 
segmentation, as shown in Figure 7. In F lg ure 7, P1 e-shaped segment S10 can be 
interpreted as being indicative of calcium, whereas pie shaped segments Sll and S12 
can be interpreted as containing no calcium. 

[0048] Other configurations of the present invention use 
segmentation of only one of the Cartesian coordinates, as illustrated in Figure 8. In 
this configuration, rectangular segments S20 and S21 are interpreted as diffenng 
iodine and/or calcium densities (concentrations) in water, assuming that no air is 
present in object 22 at the pixel of interest. (Assumptions such as this can be 
introduced in some applications of the invention by having an a priori knowledge of 
the composition of object 22.) Segment S22 can be interpreted as essentially pure 
calcium. Segment S23 is interpreted as out of range, i.e., containing no calcium, no 
iodine, and no water. 

[0049] In some configurations of the present invention, to determine 
a density of a material m, it is assumed a priori that a voxel is composed of a known 
material. More particularly, it is assumed that c m and d m are known. Then, a dual 
energy measurement is performed to obtain PmA and p mB in the relationship 
a = PmAMA + PmB p B . Using the material assumption (c m ,d m ), a new relationship is 
obtained by the relationships c m p = PmA and d m p = PmB , where p is the density of 
material m. A ratio R = PmA I PmB = cJd M is then determined, which is independent 
of p. 

[0050] For a particular pair of basis materials A and B, such as water 
and iodine, the combination of c m and d m can be assumed to be unique for any other 
material, and the ratio R can also be assumed to be unique. The ratio R is then used 
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to identify the matenal m in a particular voxel or pixel element. For example, the 
ratio R falls into one of a plurality of predetermined ranges, each corresponding to a 
limited number of target matenals. Ranges of this ratio correspond to angular 
segments of a rectangular plot, such as segments S10, Sll, and S12 shown in Figure 
6. However, in some configurations, there can usefully be as few as two segments or 
as many segments as can be accommodated within the available precision and noise 
limitations of the CT imaging system. The ranges of R are determined in accordance 
with ratios that would be obtained for various materials of interest (e.g., calcum, 
water, iodine, etc.), and the rectangular plot is divided into a limited number of 
sectors. The segments are not necessarily of equal size, but each segment represents a 
limited, predetermined range of material compositions. 

[0051] Thus, some configurations of the present invention selectively 
display only image pixels within a certain range R^ < R < ^ to display an image 
showing the locations of a selected material within an object 22. R^ and R^ are 
determined in advance for any particular material to be displayed. Other materials 
have R values outside this range and are not shown in the image. 

[0052] For example, a phantom containing gadolinium (Gd), calcium 
(Ca) and iodine (I) at specific locations was scanned by a conventional CT imaging 
system. The reconstructed image clearly showed the location of each element, but the 
composition of the spots indicative of these locations could not be determined from 
the image. By performing a multiple material decomposition, images showing only 
iodine locations free of contamination of other materials were produced by indicating 
only locations at which 0 < R < 1 . Calcium was shown in a separate image produced 
by indicating only locations at which 10<*<12. Gadolinium was shown in a 
separate image produced by indicating only locations at which 35 < R < 37 . 

[0053] It is not necessary that the basis materials or the segmentation 
of coordinates or of the ratio R be predetermined by a manufacturer, or that imaging 
apparatus 10 of the present invention be committed to specific basic materials and/or 
restricting the set of materials that can be identified. Instead, in some configurations 
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of .he present invention and as noted above, a nser is ab)e to select his or her own set 
of basis materials and generate a composition matrix as described in connecuon wdh 
Figure 4 A.so, althongh preprogrammed sections and/or ratio ranges can be 
provided ,o a user for one or more sets of basis materia!*, some configurations of the 
present invention permit a user to specify his o, her own segmentation. The user can 
determine such segmentations empiricaf.y or otherwise, to suit his or her own need, 
ta cases in which the user chooses the basis materials, ,h. basis materia* are referred 
,o herein as being -preselected" if the selection is done prior to the scanning of an 
object On the other hand, the term "postselecteti" refers to a selection of basis 
materials made after the scanning of an object. The tern, "selected" is used herem to 
broadly encompass either case. The determination of segmentation may also occur 
before or after a scan is performed. When broadly referring to either case, the tetm 
■•determine*' is used, rather than the less inclusive "predetermmed or 
"postdetermined." 

[0054] The usefulness of various configurations of the present 
invention is revealed in Figures 9 and 10, which show true lumen 500 shape and 
calcium blockage 502 in a simulation phantom obtained by multiple material 
decomposition over a complete range of Ca/I ratios. By contrast, and re emitg to 
Figure 11, conventional CT imaging is essentially unab.e to show simulated calcium 
blockages in lumens 500 of the simulation phantom. 

[0055] In still other configurations of the present invention, a gray 
scale image is useftrfty displayed in accordance with values of 
P = (^'c,+P.»'''.) /2 ' 

[0056] It will be recognized that configurations of the present 
invention are capable of providing quantitative materia, information. In 
configurations of tire present invention, only one material (subject to an overall system 
materia! separate sensitivity) is reconstructed with the correct density for each image, 
.hereby improving targe, material visibility and quantification accuracy. 
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Configurations of the present invention also can prov.de decomposition into more 
than one material. 

[0057] The present invention is useful for various applications that 

decomposition, tissue segmentation, or material quantified 

[00 58] While the invention has been described in terms of various 
speclfl c embodiments, those skilled in the art will recogmze that the invenuon can be 
practiced with modification within the spirit and scope of the clavms. 
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